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FOREWORD 


A part of the general philosophy associated with the investigation reported 
herein was the transfer of the NASA plume chemistry technology to the U.S. Navy 
to establish an in-house capability in the Navy so that any required subsequent 
analysis of similar rocket motors can be performed. The procedures, computer 
programs, users’ manuals, and input are at the Naval Surface Weapons Laboratory , 
Dahlgren, Virginia. 

This study was accomplished in response to a request from the Naval 
Surface Weapons Laboratory as given in Navy Work Order Number N<;0!)21-7(;-MP- 
H11US3. This report has been reviewed by J. W. Hill of the Naval Surface Weapons 
Laboratory and approved for publication. It contains no classified information. 
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THERMODYNAMIC AND CHEMICAL PARAMETERS 
OF THE EXHAUST EFFLUENTS FROM THE 
HARPOON BOOSTER MOTOR 

I. INTRODUCTION 


Analytical techniques were utilized to predict the thermodynamic and 
thermochemical interactions of the HARPOON rocket exhaust plume with the 
ambient atmosphere to ascertain the resulting chemical species and their 
respective concentrations so that the localized influence on the air quality could 
be assessed. 

When a solid propellant missile such as the HARPOON is launched, the 
hot exhaust effluents form a cloud behind the vehicle. The plume and the ambient 
atmosphere interact in a complex manner that is dependent in part on vehicle 
parameters such as chamber pressure, propellant ingredients, trajectory, and 
the prevailing ambient atmospheric conditions [1J. Without analytical predic- 
tions, it could be postulated that the elements in the propellant ingredients and 
the air form a voluminous number of ehemical species — some of which are 
highly toxic even in small concentrations. Thus, an analytical prediction of the 
thermodynamic and thermochemical interactions that occur in the exhaust plume 
is utilized to eliminate those species that cannot be formed. The analytical 
prediction is obtained utilizing the following analyses techniques: 

a. Thermodynamic Analysis (Step 1) — This analytical technique is used 
to examine the exhaust constituents to determine if temperatures and pressures 
are conducive to promote the formation of various species. 

b. Finite-Rate Chemical Analysis (Step 2) Using the candidate 
species selected and their postulated chemical reaction#, this analysis technique 
examines the kinetics of the flow field to determine if the residence time of the 
constituents is long enough for the postulated reaction to occur. Thus, the 
analytical predictions provide the species present and their concentrations. 

The results are applicable only for the rocket motor under investigation 
since changes in propellant ingredients, chamber temperature and pressure, 
or in the internal flow field will result in different exhaust products. 
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The analytical techniques discussed in this report serve two purposes. 
First, the exhaust species that are present and their concentrations in the plume 
are identified. If all the identified species are nontoxic, then no further action 
is required. However, if potentially toxic species do exist, it may be desirable 
to monitor the exhaust effluents. A secondary purpose of the analysis is to 
identify species that are easily monitored and determine their constancy in 
relationship to other species that are of primary interest. A hypothetical 
example is that both sulfur and chlorine compounds are present and the sulfur 
compound is easily and accurately measured, whereas the chlorine compound is 
difficult to measure. If the analysis indicates that the ratio of the sulfur com- 
pound concentration to the chlorine compound concentration is fairly constant, 
it is sufficient to measure only the sulfur compound concentration. The use of 
the naturally occurring tracer will allow accurate measurements to be made 
easily and inexpensively. 

The analytical techniques discussed in this report are applicable not only 
to the HARPOON booster but to any solid propellant rocket-powered missile. In 
addition, the computer codes, the necessary users’ manuals, and the input and 
output of all of the computer runs generated during this investigation have been 
delivered to the appropriate personnel at the Naval Surface Weapons Center, 
Dahlgren, Virginia. 

A general discussion of the overall analysis scheme that was utilized in 
the investigation is given in Secti' n II. Section III is a discussion of the thermo- 
chemical analysis and the e^rcr in the thermodynamic data utilized in portions 
of the analysis. Section IV is a discussion of the nozzle analysis and after- 
burning analysis performed. Section V gives the results, conclusions, and 
recommendations. 


II. GENERAL TECHNICAL APPROACH 


The general approach used in this study was as follows. First, a one- 
dimensional thermochemical analysis of the propellant was performed. Second, 
utilizing these results for the species, pressures, one-dimensional combustion 
product velocities, etc., a more realistic two-dimensional, two-phase analysis 
of the HARPOON motor nozzle was performed for the exit plane value of species, 
temperatures, pressures, and velocity. Finally, utilizing the results from the 
two-dimensional analysis, the one-dimensional thermochemical analysis was 
redone to obtain a realistic estimate of the species present at the motor nozzle 
exit and their concentration. Simultaneously, a study was performed to determine 
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the accuracy of the thermodynamic data utilized In the study, which is primarily 
obtained from spectroscopic constants. An estimate was made of the effect of 
the thermodynamic data accuracy on the species concentration. Using an 
available two-dimensional, finite-rate chemical mixing analysis, the plume 
emanating from the HARPOON motor nozzle was analyzed to obtain the species 
and their concentration and the local static temperature as a function of axial 
and radial distance from the nozzle exit. For some of the trace chemical 
species for which there were no data available on possible finite-rate reactions, 
the concentration of the trace species was obtained by assuming an infinite rate 
for the reaction, i.e., an equilibrium calculation was made which yielded the 
largest possible concentration that could be present. 


III. THERMOCHEM I CAL ANALYS I S FOR EX I T PLANE 
SPECIES AND CONCENTRATIONS 

A. Background 

The initial step in determining the HARPOON booster exhaust products 
that are to be employed in the air quality assessment is to eliminate as many 
candidate species as possible by a thermochemical analysis assuming chemical 
equilibrium. Basically, in this analysis the chamber and nozzle temperatures 
and pressures are used to determine the chemical reactions which can oc r ur 
without considering reaction times, hence obtaining the maximum number of 
species and their concentrations that can be present in the exhaust effluents. 

To explain in greater detail, a summary discussion of the One-Dimensional 
Equilibrium (ODE) code will be given in terms of how it functions and its data 
requirements. In this section, specific attention will be given to the effect of 
the dynamics of the rocket or. the species concentration and the treatment of 
competing reactions. The reliability of the spectroscopic data used in ODE will 
also be discussed. 

To determine what species will be present and in what concentration, a 
computer simulation of rocket firings was conducted with the one -dimensional 
thermochemical nozzle analysis program (ODE). The input for the analysis 
consists of the solid propellant composition and heat of formation, the HARPOON 
motor geometric specifications, and the chamber pressure history. The ODE 
computer program [2] determines the chemical species and their concentration 
at chemical equilibrium for a given motor condition, i.e. , chamber pressure 
and- area ratio. The resulting species of the equilibrium calculations were used 
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as input to tlio nozzle analysis and afterburning calculations. Inherent in the 
calculation of the chemical species and their concentration are the potential 
errors in the thermodynamic data which produce some uncertainty in the 
concentration of chemical species. 

•* 

The chemical species present in the combustion chamber and their con- 
centrations depend on the combustion temperature. The combustion temperature 
depends on the enthalpy of the propellant ingredients and the chemical specie ' 
formed in the combustion process. The solution of such a cyclic interactio’ ; - 
calculated by a one-dimensional thermochemical program. One of the requuod 
inputs to such r program is the propellant ingredients and their enthalpy. Ihe 
HARPOON booster propellant formulation and enthalpy used is given in Table i. 
It should be noted that the propellant contains 78.5 percent of a highly reactive 
compound, ammonium perchlorate, as an oxidizer. 1 The use of such an 
oxidizer implies that the combustion products should contain significant 
quantities of chlorine species. A motor parameter which is significant when 
calculating the exhaust species is the expansion or area ratio of the nozzle. 

Table 2 [3] gives some of the important geometric parameters of the HARPOON 
booster. An additional discussion of the significance of some of the parameters 
as they apply to the finite-rate analysis is given in Section IV. 

TABLE 1. HARPOON PROPELLANT COMPOSITION 


Ingredient 

Weight Percent 

Ammonium Perchlorate 

78.5 

Alumi nu in 

G.O 

Copper Chromite 

2.0 

Sulfur (Flowers) 

0.1 

Binder (PU) 

13.4 

Enthalpy of propellant (cal /mode) 

-27i 


The chemical species computed in the one-dimensional equilibrium 
thermochemistry program are based on the chemical composition of the pro- 
pellant and the chemical compounds in the program data base. First, all 
species that are chemically possible for the propellant under consideration are 
identified. Based on thermochemical data, the compounds listed in Table 3 [ 1| 

1. Private communication from Dr. Glen Moore, Naval Surface Weapons Center, 
Julj 1975, 
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TABLE 2 . HARPOON MOTOR PARAMETERS 


Throat Diameter 

2,23 in. 

Throat Area 

3.91 in . 2 

Nozzle Diameter 

6 ,‘H' in. 

Nozzle Area 

32.17 in . 2 

Nozzle/Area Ratio 

8.24 


TABLE 3. CHEMICAL SPECIES CONSIDERED 


Cu 

NH 

H 

COS 

ai 2 o 2 

A1 

S 

HCO 

C 2 C1 2 

C 

AlClj(L) 

so 2 

h 2 o 

C 2 0 

CH 

AlO 

Cr 3 C 2 (S) 

NO 

CIO 

CH 4 

A10 2 K 

CuCl(L) 

n 2 

CR(S) 

CO 

A1 2 0 3 I s) 

Cu(L) 

0 

HALO 

co 2 

CC1 

A1C1 

SH 

hno 2 

C 2 I1 

ch 2 

AICI 3 

S0 3 

I 1 2 0 2 

c 3 0 2 

CN 

A10C1 

Cr 7 C 3 (S) 

NOC1 

C10 2 

COC1 

A1 2 C1 6 

CuCl 

n 2 o 

CR(L) 

cs 

ai 2 o 3 (l) 

CuO 

OH 

HCl 

c 2 n 

CC1 2 

A1C1 2 

SN 

1 I 0 2 

Cl 

ch 2 o 

A1N(S) 

S0 2 C1 2 

1I 2 S 

Cl 2 

CNN 

AlOH 

CrCl 2 (S) 

CM 

O 

>5 

Cr 

COCl 2 

ai 2 o 

CuH 

nh 2 

HCN 

cs 2 

C(S) 

Cu 2 

0 2 

h 2 

c 2 n 2 

CC1 3 

A1C1 3 (S) 

SO 

NHj 

C1CN 

CH 3 

AIN 

s 2 

no 2 ci 

C1 2 0 

cn 2 

aio 2 

Cr 2 0 3 (S) 
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were selected as the basic candidate set of species. The analyzer then deter- 
mined which of the various chemical species existed in the simulated HARPOON 
motor. The thermochemical data required for this analysis included the heat 
capacity (C ), entropy (S), enthalpy (H), Gibbs free energy (G),*and equilib- 
P .at 

rium constants (K ) as a function of temperature. Finally, these data*were 
P 

used by the thermochemistry program to compute the concentrations of the 
chemical species in the chamber, throat, and exit plane based on the propellant 
enthalpy and the temperature and pressure that occurred locally. , 

As previously noted, the concentration of a speqjes is a function of the 
chamber pressure and temperature. To demonstrate this, the first 2 sec of 
data for three trajectories of the HA7 ’OON booster for each oil three HARPOON 
launcher configurations, the Canister, Tarter, andASROC, are tabulated in 
Table 4 [5J. The three different launcher configurations provide different angles 
of launch and initial launch positions relative to the ship% deck. This gives 
different angles of the missile centerline and ranges relative to the desk as a 
function of time. A chamber pressure history [ 6 J is also provided ir* this table. 
Both sets of data were analyzed a function of time. The analysis yielded 
range, velocity, and chamber p Insure tabulations as a function of time .(Table 4) . 
Figure 1 show's the missile trajectory :or the Canister launcher. 

A series of computer runs was made to determine the c®ncentration of 
chemical species at different time periods in the initial portion of the flight of 
the missile. The data for several pressures are given in Table 5. The species 
concentrations are at the 1-D exit plane. The 1800 psia is indicative of the 
chamber pressure just after the first motion out of the launcher. Table 5 lists 
only those compounds which contain significant concentrations (above 10 ppm). 

It can be seen that most species concentrations changed very little from the 
pressure at the initial phase to peak chamber pressures 1 } Ammoma (NH 3 ) , 
carbonyl chloride (COC1) , and hydrogen cyanide (HNC) increased approximately 
20 percent in concentration as the pressure increased to the peak. This is 3 
small change over the 300 psia pressure differential. It is understandable 

t cause the equilibrium constant, K , for formation of products is very 

P 

temperature dependent. Consequently’, it is not critical that one investigate the 
impact of various chamber pressures, but it is more important to examine the 
portions of the plume where the temperature increases significantly. 

In a combustion process, normally there are competing reactions between 
some of the species. The production of particular chemical species and their 
concentration under equilibrium combustion conditions at a given temperature 


r, 
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TABLE 4. RANGE, VELOCITY, AND PRESSURE VERSUS 
TIME SINCE MISSILE FIRST MOTION 
OUT OF LAUNCHER 




Time 

(sec) 

Chamber 

Pressure 

(psia) 

Range from Deck (ft) 

Velocity (ft/sec) 

Tartar 

ASROC 

Canister 

Tartar 

ASROC 

Canister 

0.1 

1800 

10.8 

19.7 

17.9 

44.9 

108.6 

115.27 

0.2 

1870 

15.4 

30.8 

30.3 

71.4 

144.0 

149.7 

0.3 

1875 

23.1 

45.9 

46.3 

99.8 

177.0 

180.9 

o;4 

1930 

34.0 

64.8 

65.8 

128.5 

211.3 

214.1 

0.5 

2080 

48.2 

87.3 

88.6 

162.8 

246.2 

246.6 

0.6 

2100 

65.8 

113.4 

114.7 

186.5 

280.3 

278.3 

0.7 

2110 

85.7 

142.9 

143.9 

226.4 

314.1 

309.5 

0.8 

2100 

110.1 

175.8 

176.3 

265.8 

347.6 

340.6 

0.9 

2000 

138.2 

212.1 

211.8 

288.1 

379.4 

371.6 

1.0 

1960 

168.7 

251.5 

250.4 

322.4 

412.6 

402.5 


1955 

202.4 

294.4 

292.1 

352.7 

445.6 

433.0 


1950 

239.2 

340.5 

336.8 

384.2 

477.7 

463.0 


1950 

279.2 

389.8 

384.5 

414.9 

510.0 

492.7 


1950 

322.2 

442.4 

435.1 

445.1 

541.4 

521.7 

^9 

1950 

368.3 

497.9 

488.6 

475.5 

570.7 

550.2 

1.6 

1945 

417.2 

556.5 

544.9 

505.2 

602.1 

578.1 


1945 

469.2 

618.1 

604.0 

534.9 

630.8 

605.2 


1940 

524.1 

682.5 

665.8 

562.7 

658.7 

631.4 

1.9 

1910 

581.6 

749.7 

730. 1 

589.8 

686.0 

656.3 

2.0 

1890 

642.1 

819.5 

796.9 


I 
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CANISTER LAUNCHER 



Figure 1. HARPOON booster trajectory. 


and pressure depend on the chemical composition of the propellant and the Gibbs 
free energy of the chemical species. To illustrate this, consider the following 
arbitrary reaction, j, with chemical species A, B, C, and D; their concentration 
will be determined by the equilibrium constant K. , where 


K. 

J 


‘ P C |C 1P D |d 
IP A ia ' P B |b 


( 1 ) 


For reaction j, 


aA + bB * cC + dD 


( 2 ) 
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TABLE 5. RESULTS FROM THERMOCHEMICAL ANALYSIS (ONLY 
CONCENTRATIONS FOR GASEOUS SPECIES) 


> 4 . 




Concentration** *1 1 

-D Exit PtiU* (ppm) 





Pr***ur# (|*U) - 

- Temperature (K) 



Sp«ci** 

1800-322* 

1870-3232 

1930-3234 

1960-3235 

2000-3236 

2080-3238 

2110-3239 

am 

6.0475+2 

6.8093+2 

5.6194+2 

5, 5288+2 

5.4126+2 

5,1936+2 

5, 1156+2 

E9 

1.3782+0 

1.3324+0 

1.2780+0 

1,2568+0 

1.2297+0 

1,1787+0 

1.1604+0 

AlCl 

1.2570-2 

1.2289-2 

1.2061-2 

1.1950-2 

1.1809-2 

1. 1537-2 

1.1438-2 

AlCl, 

1.1298-1 

1. 1268-1 

1.1244-1 

1,1232-1 

1.1217-1 

1.1187-1 

1.1175-1 

Aid) 

1.0370+0 

1,0563+0 

1.0720+0 

1.0797+0 

1.0900+0 

1,1101+0 

1.1175+0 

A10C1 

8.4671-2 

8.2805-2 

8,1293-2 

8.0561-2 

7,9621-2 

7.7819-2 

f. 7163-2 

AlOjH 

2.2437-2 

2. 1935-2 

2. 1529-2 

2, 1332-2 

2.1080-2 

2,0596-2 

2.0420-2 

CN 

4 . 174 -7 

4,243-7 

4.301-7 

4.329-7 

4,367-7 

4.441-7 

4,468-7 

co 0 
COCl 

1.0740+5 

1.0745+5 

1.0743+5 

1.0743+5 

1,0749+5 

l,0?47+5 

1,0740+5 

3.3739-2 

3.4341-2 

3,4847 -2 

3. 5096-2 

3,5427-2 

3.6076-2 

3.6315-2 

COCl, 

CO, 

1.193«^k 

1.000 

1.238-6 

1.0009+5 

1.276-3 

1.0010+5 

1,296-5 
1. 0011+5 

1.321-5 

1.0012+5 

1.373-5 

1.0013*5 

1.392-5 

1.0014+5 

Cl 

2.1589+2 

2. 1150+2 

2.0794+2 

2.0622+2 

2.0399+2 

1,9973+2 

1.9819+2 

C1CN 

8.979-0 

9.315-6 

9.603-6 

9.747-6 

9.939-6 

1.032-5 

1.047-5 

Cl, 

4.7017-1 

4.6951-1 

4.6899-1 

4.6872-1 

4.6840-1 

4.6775-1 

4.6749-1 

HC1 

1,8237+3 

1.8238+5 

1.8239+5 

1.8239+5 

1.8239+5 

1.8240+5 

1.8240+5 

HCN 

1.0504-2 

1.0905-2 

1.1249-2 

1.1421-2 

1.1650-2 

1,2108-2 

1.2280-2 


2.6441*3 

2.7126+2 

2.7698+2 

2.7981*2 

2.8351*2 

2.9078+2 

2.9346+2 

NH, 

5.2735-1 

5.4782-1 

5.6536-1 

5.7413-1 

5. 8582-1 

6,0920-1 

6.1798-1 

NO 

1.3875+0 

1,3581*0 

1.3342*0 

1.3226*0 

1.3078+0 

1.2792+0 

1 * 2689+0 

NO, 

2.771-6 

2.708-6 

2.658-6 

2.633-6 

2,601-6 

2. 541-6 

2.519-6 

SO 

5.4139+1 

5.3373+1 

5.2735+1 

5. 2421*1 

5.2009+1 

5.1204+1 

5.0907+1 

SO, 

4.7999+2 

4,7336+2 

4,6787+2 

4.6518+2 

4.6164+2 

4.5471*2 

4,5216+2 

so, 

1.7469-3 

1,7210-3 

1.6994-3 

1.6888-3 

1.6750-3 

1.6479-3 

1.6379-3 

CuCl 

1.9606+1 

1.9270+1 

1.8948+1 

1.8792+1 

1.8591*1 

1.8207+1 

1,8067+1 


t. 6.0475+2 « 6.0475E+2 


The values in the brackets are the partial pressures for the species. The a, b, 
c, and d in equation (2) are the stoichiometric coefficients or the number of 
molecules of each species. 


The equilibrium constant is related to the Gibbs free energy of the 


standard state for each chemical species by equation (3): 


K 


j 


-ag At 


» 


e 


( 3 ) 
















where 


AG « cG c + dG D - aG A - bG fi . (4) 


The G is the Gibbs free energy at the temperature in question for the standard 
state of the chemical species. 

In the actual combustion chamber many reactions occur concurrently. 
Some are competing reactions which require some of the same species as a 
reactant. The thermochemistry program simultaneously calculates the con- 
centrations of the many possible species by minimi ziilg the Gibbs free energy 
for the combustion chamber under equilibrium conditions. Thus, the chemical 
species computed in the one-dimensional equilibrium thermochemistry program 
depend on the chemical composition of the propellant and the chemical compounds 
in the program data base. The thermochemical data for compounds listed in 
Table 3 were the basic set of data which determined the various chemical species 
that can exist in the simulated HARPOON motor. 

In actual practice, one must consider all potentially applicable chemical 
species. However, because the computer program is limited to handling 100 
species, those species with molecular fractions less than 10 ~ 18 are deleted. 
Nevertheless, one must always prove that the deletion of a species had no sig- 
nificant effect on the concentrations of the other species in the system. 

A small study was conducted by omitting species that had negligible 
concentration in a run and reinserting species that had been previously omitted. 
No species concentration changed, thus proving the validity of the technique. 

The concentration of the species depends on the equilibrium constant K, 
which is dependent on the temperature via equation ( 5) : 


d In K _ H 
dT " kF 


( 5 ) 


where H is the enthalpy of the reaction being considered [7] . Consequently, 
it is not critical that one investigate the impact of various pressure chambers. 
It is noted from Table 5 that the temperature increase is 10 K over the entire 
chamber pressure range. This small increase in temperature causes a 
corresponding small increase in concentration of a few species. 
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In Paragraphs m,B and m.C, an error analysis of the resulting species 
and their concentrations will be presented. However, before the error analysis 
is discussed, the source of the thermodynamic data used together with its 
accuracy should be explained. The thermodynamic data needed to calculate file 
gaseous chemical species in the plume are readily computed via partition functions 
from spectroscopic constants. The thermodynamic properties of liquids and 
solids must be obtained from experimental measurements. The thermodynamic 
properties are normally expressed as functions of the internal partition function 
Q, i. e. , 


C P = T 2 d 2 ln Q 2Td(ln Q) 5 

R dT 2 dT 2 


where the superscript o is standard state of the heat capacity at a constant 
pressure C , R is the universal gas constant, and T is the temperature. Then, 


S T = Td(lnQ) 
R dT 


h- InQ + fin M + flnT + S + 

M Z C 


5 

2 


CO 


where S is the entropy and M is the number of molecules, 


H T ~ H 0 _ TdlnQ 5 

’ “ RT = dT 2 


( 8 ) 


where H^, is the enthalpy at a temperature, and 





a e 

In Q + — In M + — InT + S 
z z c 


( 9 ) 


where G^, 1 b Gibbs free energy at a temperature and 
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The internal partition function Q is represented as 




m=l 


( 10 ) 


where Q m is the internal partition function for the mth electronic state and L is 
the number of electronic states. 

The molecular partition function Q m is described in equation (11): 


Qm _ m ^m _m _m _m m _m 

= Q Q Q Q Q Q Q 

^e v ^0 x 


( 11 ) 


where Q g , Q y , and are the electronic, vibrational (harmonic oscillator) , 

and classical rotational contributions. The other quantities are as follows: 

rotational stretching, Q*P; low temperature rigid rotation, Q™; Fermi resonance, 
m ^ m 0 

Q ; and both anharmonic and rotation interaction, Q , are corrections. The 
a; c 

last four terms are correction terms for various interactions not considered 
classically. Partition functions are classically thought of as a way to partition 
the internal modes of energy into rotation, vibration, and electronic modes. 
These energy modes are coupled and depend on each other. 


The partition functions themselves can be calculated very accurately for 
a given electronic, vibrational, and rotational state, but molecules of particular 
temperatures will contain many electronic, vibrational, and rotational states. 
Thus, techniques for terminating the terms used in the calculation of the 
thermodynamic properties have been developed [8,9] . 
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B. Error Analysis of Atomic Species 

The uncertainty of the thermodynamic properties of the atomic species 
is now considered by examining the spectroscopic equation governing these 
thermodynamic characteristics. 

In a monatomic molecule, the electronic energies are the only internal 
energy contributions; therefore, equation (10) applies only to the: electronic 

o 

contributions. The use of the electronic states to compute the entropy (S ) , 
heat capacity (C^), and enthalpy (H) for the atomic species is limited by deter- 
mining the number of electronic states to be considered. At the temperature 
considered in the HARPOON combustion chamber (<3000 K) , the thermodynamic 
properties are accurately described, because many upper electronic states are 
not heavily populated at these temperatures. In some atoms, such as Cr, there 
are many electronic states to consider, although, again, at these temperatures 
this does not present a problem. The excited electronic states are also fairly 
well known for atoms. An error of 10 percent is considered an upper limit for 
the thermodynamic properties for any atomic species. 

m 

C. Error Analysis of Diatomic and Polyatomic Species 

The spectroscopic constraints for the vibrational, rotational, and correc- 
tion terms are not as well known as for electronic states. This is particularly 
true when these properties must be determined at many different excited elec- 
tronic states for an infinite list of compounds. In all cases the correction terms 
are included in the Joint-Army -Navy-NASA -Air Force (JANNAF) tables. If the 
spectroscopic constants were not known for specific compounds, they were esti- 
mated from analogous works with similar compounds. There are also some errors 
in assignment of states by different researchers. The errors for these com- 
pounds are listed in Table 6 for entropy, enthalpy, and spectroscopic constants; 
an overall error estimate is also given. 


D. Summary of Error Analysis 

The analysis of errors in the thermodynamic data is straightforward in 
many cases but somewhat subjective in others. In general, the overall error 
approximation reflects the highest error for a particular thermodynamic 
parameter. Since the parameters are all Interrelated, errors tend to propagate 
through all of the parameters. It should be noted that these are probably worst- 
case error bounds. The error in NH 2 is the most severe. The values which are 


13 



TABLE 6. ERROH ANALYSIS 


Species 

s° 

(%) 

AH f ° at 

298 K 
(%) 

Dissociation 

Energy, 

D 

0 

(%) 

Vibrational 

Levels, 

v i 

(%) 

Appropriate 
Overall 
Error (%) 

Comments 1 

A1C1 

■ 

8 

1 


8 


AlCi, 

III 

7 



7 


A1C1 3 


<1 

■■ 


1 


I AlO 


12 



12 


AlOj 

3.4 

11.3 

■■ 


10 


AlOjH 

2.5 




5 


AlOH 

1.0 

7 

■H 


7 


AljO 

2.4 

16 


15.8 V x 

16 


AljOj 

6.0 

9.6 

■ 


9.6 


AljOj(S) 

AljOj(L) 

<1 

<1 



1 

Very well studied 

C CIO 


66 



66 


Cl 





5 

Molecule studied in 
great detail but many 
higt’t v electronic 
states 

Cl 2 





3 

Some errors in 
electronic states 

C1 2 0 


2.8 



2.8 


Cr 

<1 

1 

il 



5 

Many electronic 
states to consider 

Cu(L) 


1 



1 

Very well studied 

Cu 


<1 



10 

Many electronic 
states and heat of 
sublimation 10 
i percent off 

CU| 


<1 


! 

i 

10 

i 

CuO 

<1 

13.8 



13.8 

Poor electronic 
assignments 

H 

ITN 


4 



1 

4 

Very well 

SH 


3.6 



5 

Analogous molecular 
parameters used [ 

h 2 





1 

Very well studied 

h 2 n 


3.2 



30 

Higher values 

could be off by factor 
of three 

h 2 o 





1 

Very well studied 
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TABLE 6. (Concluded) 





Dissociation 

Vibrational 





AH f ° at 

Energy, 

n 

Levels 

v 

Appropriate 



s° 

298 K 

0 

v i 

Overall 


Species 

(%) 

(%) 

<%) 

(%> 

Error (%) 

Comments 

HA 





1 

Very well studied 

HjS 


4.8 



5 

Well Studied 

NHj 


1 



1 

Very good data 

NO 


<1 



1 

Very few electronic 
states 

NO, 


%6 


m 

2.5 






1 

Very good data 

N,0 


<1 



1 

Very good data 

O 





1 

Very good data 

O, 

<1 




1 

Very good data 

so 

<1 

25 



20 

All other data 


Vi 





fairly good 

SO, 

<1 

<1 



1 

Good data 

so, 



I 


1 


s 

<1 

<1 



1 


S* 


<1 



1 


OH 


3 



3 

Some errors below 
400 K 

CO, 

<1 

<1 



1 


COS 


<1 



1 

Accurate spectro- 
scopic data 

CN 

<1 

2.4 



2.4 

Accurate spectro- 
scopic data 

NOC1 

<1 

<1 



1 


HC1 


<1 



1 

Accurate spectro- 
scopic data 

CuCl 


1.8 



1.8 


ch 4 

<1 

<1 



1 


CH, 

<1 

<1 



1 

1 


CH,0 

<1 

5.4 


m 5.4 

Accurate spectro- 







scopic data 

CH, 

1 

1 



1 


HCO 

<1 

1.9 



1.9 


CO 

<1 

<1 



1 

Very well studied 

HCN 

<1 

6.0 



6 

Fairly good data 

AlOCl 


6.0 



6 

Fairly good data 




•» 
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calculated using spectroscopic data are, in general, accurate. For liquids and 
solids, the data can be measured only experimentally, and usually larger errors 
are found. It should be noted that it is generally agreed that the calculation of 
thermodynamic properties via partition functions accurate or more 
accurate than allowable experimental errors. Generally, the errors are small 
and should give small errors in concentrations. Errors in concentration of 
100 percent should be considered extreme cases and then only as a margin of 
safety for potentially hazardous species at high dosagjes. This condition is true 
because of the relationship 

I# 

-A G t 

2.302 RT = 1Og 10 K p * 

where AG^, is the Gibbs free energy difference between reactants and products 

and Kp is the equilibrium constant of the chemical reactions at equilibrium. 

^ is the ratio of the concentrations of the reactants and products. The log 

term magnifies any error in AG^,, but all Gj must be in error in opposing 

directions for serious errors to present complications. As a result of the many 
competing reactions, serious errors in certain species become diluted to some 
degree. The data for error analyses came from the JANNAF Thermochemical 
Tables, Volume I, II, HI, and KT and their references [10]. 

Thus, the species listed in Table 5 will be the basis of the next step in 
the analysis. 

€ 

IV. SOLID ROCKET MOTOR NOZZLE AND PLUME ANALYSIS 


The assumptions and resultant species for the nozzle and plume are 
presented in this section. The results obtained in Section III are used as the 
inputs to this analysis. 

The assumption and theoretical analysis are very important in arriving 
at the species and their concentrations in the plume for the following reasons. 
First, in the chamber analysis the system was well defined in that all the 
constituents and the fluid dynamics involved in the reaction were well defined. 
However, the constituents and the fluid dynamics of the rocket exhaust plume 
are not as easily defined. More specifically, air is now entrained into the 
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plume and shock fronts are set up. Second, very little experimental data on the 
species have been directly collected »nd analyzed from a solid rocket motor 
plume to support the theoretical analysis* This lack of empirical data for the 
hot exhaust plum»exists due to the complexities of obtaining meaningful measure- 
ments from the plume. Therefore, the only way to verify the accuracy of the 
analyses techniques utilized in this study is that the performance of a solid 
rocket motor, whichdepends oft the species present and their concentration* 
can be predicted to a high degree of confidence. Secondary techniques have been 
used (s*g. , sampling the stabilized exhaust cloud) which afford indirect verifica- 
tion of the analysis; however, this sampling provides only a ooarse verification 
at best. 


A. Background of Analysis Logic 

9 

The order of our analytical procedure for a solid rocket motor was to 
first perform a one-dimensional thermochemical analysis and then use the 
results of that analysis as the basis for a two-dimensional nozzle thermo- 
chemical and aerodynamic analysis. For solid propellants with aluminum in 
the formulation, a two-dimensional, two-phase nozzle analysis is required. 

Since our primary interest is in the species In the plume (and since good plume 
experiments do not exist) , we must use indirect verification from a solid rocket 
motor performance analysis for an indication of species concentration accaracy. 

{3t 

Figure 2 shoes the schematic used for the motor, nozzle, and plume 

analysis. The first portion of the figure shows the solid propellant motor 

burning at some chamber pressure, I , aid at chamber temperature, T , 

c c 

and'some of the principal (but not necessarily all) speciee of Interest that are 
formed in the combustion process. This combustion analysis is valid up to the 
nozzle throat, where two-dimensional effects start. From the throat to the 
nozzle exit plane, a two-dimensional, two-phase nozzle analysis is utilized. 

From the nozzle lip downstream, a region exists which is primarily controlled 
by the mixing of the ambient air, and chemical reactions which are occurring 
at a finite rate are analyzed. Therefore, we have to perform a finite-rate mixing 
and afterburning analysis. 


B. Analysis of Losses 

The losses considered in the investigation are given in Tatte 7. The 
normal flfezle losses (divergence, two-phase Cow, and changes in the tempeft** 
ture and pressures in the nozzle) affect the species. The boundary layer jfoss 
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Figure Z* Solid rocket motor, nozzle, and plume 


TABLE 7, LOSSES CONSIDERED BY PERFORMANCE ANALYSIS 


• Divergence 

o 

• Two Phase 

• Boundary Layer 

• Energy Release 

• Submergence 

• Kinetic 

• Combustion in Nozzle and 
Plume 


Two dimensionality of nozzle 
Velocity and temperature lag of particles 
Effect of friction and heat transfer 
Coupled C* (subsonic) 

Nozzle entrance effects 
Chemical nonequilibrium 
Afterburning 


analysis determines the inviscid conto'ir that is analyzed; the ei. release lose 
is an indication of the combustion efficiency of the motor. The i . loss is as 
indication of the rate of reaction in the nozzle* In the case of the HARPOON booster, 
the motor is operating at high chamber pressures on the order of 18CK) to 2100 psi 
[11), and kinetics is not an important factor. When the flow' is expanded in the 
plume, the pressure and temperature drop rapidly. The kinetics of the reactions 
in the plume ore important, and a finite.-rate mixing and after-burning analysis is 
required. A loss that is not considered but which does occur in some motors is 
the loss due to nonuniform entrance conditions. Since the HARPOON grain is axi- 
symmetric and has a conventional design, this loss need not be considered in the 
particular case of the HARPOON. A variety of analysis tools is available (Table 
8). The variety of one-dimensional programs gives an indication of the losses and 
generates input to the two-dimensional programs. The results of the two-dimen- 
sional programs are then fed back to the one-dimensional programs to get species 
concentrations. These species and exit conditions are used as input to the low 
altitude plume program (LAPP) code for an afterburning analysis. 

Table 9 gives the input requirements of the two-phase nozzle analysis 
program. A chemistry model has to be chosen, physical constraints (e.g. , the 
nozzle contour and external ambient conditions) must be determined, and the 
particle size and mass flow distribution must be chosen. The other parameter 
requirements are usually set and are not changed from case to case. 

o l. Losses Due to Two-Phase Flow . An important consideration is 
where to perform a one- or a two -ph ace flow analysis. One-phase flow means 
that the rocket exhaust can be treated as containing only gaseous species, where- 
as two-phase flow means that rocket exhaust must be treated as being composed 
of both gaseous and solid species. 
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TABLE S. ANALYSIS TOOLS AVAILABLE 


• ODE 

• OD2P 

• ODN ' 

• QDL’PK 

• TDE 

• TPKNA 

• TDK 

• TB’L 
© AGP 

• SCR 

• LAP"" 


One-dimensional theoretical equilibrium performance 
One-dimensional, two-phase perfect p as nozzle analysis 
One-dimensional, one-phase real gas kinetic nozzle analysis 
One-dimensional, two-phase real gas kinct|c nozzle analysis 
Two-dj . ^ensional, one-phase real gas nozzle analysis 

j 

Two-dimensional, t\ -phase real gas nozzle a:.alys*s 
Two-dimensional, one-phase real gas kinetic nozzle analysis 
Turbulent, compressible boundary layer analysis 
Thermodata preparation 

Two-dimensional, two-phase, real gas nozzle analysis 
Plume afterburning 

: Al. 


TABLE 0. TDDD a INPUT REQUIREMENTS 


O 


Che«mistry 

• Frozen (ideal gas) 

• Equilibrium ° O 

(?k 

• Equilibrium /frozen 

• Finite rat^ ® 

Physical Constraints ® 

• Nozzle contqjur and external conditions 
Particle Information 

• Size and mass flow distributions (enthalpies, C , 

heat of formation (JANNAF) ^ 

• rifts ity,. 

Gas/Particle Interface Data 

• Drag coefficients c 

«> Heat transfer coefficients 


-o-rP- 


<0 


’a. Two-dimensional, two-phase nozzle analysis program 








Samples collected at rocket firings (Fig. 3) clearly show the presence of 
aluminum oxide 1 12 J . 

As can be seen, a number of spheres of varying size have been collected. 
The largest particle seen at lliis test was on the order of 21) pm in diameter; the 
smallest particles, which are small white dots seen alone on the figure and 




occasionally on tlu* larger particles, ;uv approximately 0. 1 pin in diameter. 
I'bese small particles (alpha aluminum oxide) an* spherical; own tlu* larger 
olios which appear hollow arc sphoiical in shape. In Figure I one can observe 
:• white tuzz aluminum oxide on top of the alpha aluminum oxide (00 pm diameter) 
This white tuzz material is an example of gumma aluminum oxide which results 
from some small solid rochet filing and probably would not be found in the large 
II \K1VH>N booster | »2|. 
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Figure I. Gamma aluminum oxide particles. 


rhe determination of tlu* size, type, anil distribution of tlu* aluminum 
oxide is ot primary importance to the analysis. The basic assumptions that ai 
ma le for particle flow are given in Table 10. The assumption of a spherical 
particle is probabh good since this is typically what is seen in our samples. 
\ll the other assumptions are standard. The analysis does not look at every 
particle. Only a simdl number of discrete particles are analy zed, and a repri 
sentative size and weight flow of each of these particle sizes is chosen. At tlu 
present time we can consider. 10 different particle sizes. 






TABLE 10. BASIC ASSUMPTIONS FOR PARTICLE FLOW 


• The particles arc spherical in shape. 

• The particle internal temperature is uniform. 

• The pas and particles exelianpe thermal energy by convection 
and radiation (optional). 

• The pas obeys the perfect pas law. 

• The pressure of the pas and the drag oi the particles contribute 
to the force acting on the control volume. 

• The pas is inviscid except for the drag it exerts on the particles. 

• There are no particle interactions. 

• The volume occupied by the particles is negligible. 

• There is no mass exchange between the phases. 

• A discrete number of particles, each of different size or chemical 
species , is chosen Li represent the actual continuous particle 
distribution. 

• The particles are inert. 


The average aluminum oxide particle diameter plotted against throat 
diameter in inches for various aluminum levels from 14 to 20 percent is shown 
in Figure 5. The majority of the data is at 10 percent aluminum, which includes 
Titan, Minuteman, ;uid Polaris type motors. From intuition and some pre- 
liminary data, the HARPOON particle size was chosen to be approximately 
12 pm diameter for a 2.23 in, throat diameter. 

2. Two-Dimensional 1 .osses . For a two-dimensional nozzle analysis, 
the nozzle must be modeled in great detail. Figure 0 shows all the various 
angles and physical parameters that are required for valid analysis. 

There are two primary techniques wliich can be used to analyze two- 
dimensional, two-phase nozzles: the mesh point technique, which was written 
by Kliegel and Nickerson at TRW Systems In the early GO's, and the streamline 
normal technique, which was developed in the mid-GO's 1 1 3 J . For this 
particular an illy sis, the mesh point construction method of Kliegel and 
Nickerson was chosen. Figure 7 shows a schematic of the calculation technique. 
The nozzle is analyzed by having a finite difference mesh which is calculated 
repetitively from the initial value line to the exit of the nozzle. Knowing the 
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Figure 6. Nozzle parameters. 
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Figure 7. Calculation technique schematic. 
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values of the flow variables and the species at two reference points, the value 
of all the parameters can be calculated at a third point which is the intersection 
of the cones of influence of these two reference points. 

The schematic of the thermochemistry table construction is shown in 
Figure 8. The table is constructed using the data given in Table 5. Since the 
two-phase nozzle flow is nonequilibrium and dissipative, the stagnation pressure 
at every point in the nozzle is different; therefore, if the species are to be 
calculated, a table which uses two calculable properties as independent variables 



PARTICLES 
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*2 >*1 


H G 2 * h g , 


RESULTS FROM 
Q*Wp 



S 


Figure 8. Schematic of thermochemistry table construction. 
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must be constructed. The present technique uses enthalpy and entropy. Knowing 
the enthalpy and entropy of any field point in the nozzle and utilizing a table 
look-up routine for the calculated stagnation pressure, the species at that static 
pressure can be determined. The species concentration determined by this 
technique may or may not be significantly different, but it has more validity than 
the value determined from a one-dimensional equilibrium calculation or a single 
variable table look-up routine. 

The results of this type performance analysis are presented in Table 11 
for seven motors (four NASA and three Navy motors) [14] . The three Navy 
motors tested were involved in the C-4 program. With the exception of the LS-2 
motor, which had difficulties because of nozzle erosion, all of the analyses of the 
performance are within ±0.5 percent. This is certainly adequate and is a good 
indication that the species calculations are valid. 

TABLE 11. PERFORMANCE ANALYSIS RESULTS 


Motor Parameters 

F'V-4 

Jupiter 

Orblter 

JPI.-ATS 

Old 

Surveyor 

LS-1 

LS-2 

DEMO-1 

Geometry 








c 

47.4 

40 

35 

46.22 

14.8 

14.7 

10.56 

R t - in. 

1.16 

0.7245 

2.04 

1.765 

2. 07 

2.07 

2.02 

Rc Rt 

2.0 

1.5 

3.0 

2.1 

4.0 

4.0 

4.0 

0 - deg 

20.0 

30.0 

26.0 

34.0 

25 

25 

26 

max 








'’exit - de * 

20.0 

11. G 

10.0 

20.49 

19 

19 

23 

&e — deg 

0.0 

18.4 

16.0 

13.51 

6 

6 

3 

L - In. 

19. 17 

15. 57 

31.83 

21.18 

15.8 

15.8 

11.85 

Ballistics 








P c — psia 

704 

143 

200 

Oil 

778 

950 

1003 

T — °R ! 

c 

6017 

5827 

5654 

6008 

6305 

6305 

6305 

C*- ft sec 

51 GO 

5164 

5097 

5195 

5289 

5288 

5288 

1 

I at c — sec 

vac 

315.3 

324. G 

306.4 

319.9 

301.38 

301.8 

293.21 

w — lb sec 

19.94 

1.55 

17.06 

31.80 

66.24 

80.7 

80.64 

p ! 








Performance 








! 

A Is 2D-2P-sec | 

20,95 

37.65* 

15.96 

25. 89 

23.92 

22.4 

23.97 

AI« ER — see 

3.75 

0. /4 

3.45 

10. 24 

0 

0 

0 

i 

A Is KIN — sec 

0.0 

1.95 

1.00 

0.0 

0.0 

0,0 

0.0 

AIs BL — sec 

3.15 

1.65 

5.95 

3.52 

2. 1 

2.1 

1.97 

‘vac pred “ ® ec 

297.5 

274.2 

280. Cl 

280.3 

275. 36 

277.3 

u 

267.37 

I delivered — sec 

vac i 

2 HO. 4 

27b. 0 

280. 2 

279.3 

275.2 

274. 3 b 

267.7 

‘vac preC — ‘vac del 

1.1 

-1.9 

-0.16 

1.0 

♦0.1 

-3.0 

-0, 33 

I deviation 

0.35 

-0.55 

-0.05 

0.31 

+ 0.04 

-1.0 

0.12 


a. Loss with single particle size; calculation with particle size distribution showed shocks in nozzle. 


b. Silica phenolic evit cone — severe erosion. 
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As further evidence of the validity of the technique, Figure 9 shows the 
ratio of chamber pressure to wall pressure versus area ratio for a 16 percent 
aluminum solid propellant fired in two motors at NASA/MSFC. The data points 
and error bands are shown as a circle with a band. The analysis shows a good 
comparison with the data. Another comparative pressure curve is shown in 
Figure 10 with the same motor fired with a 2 percent aluminum solid propellant. 
Again, good results were obtained using the analysis technique. Thus, since the 
HARPOON motor has 6 percent aluminum in the propellant and the analysis has 
bounded the problem at both the 2 and 16 percent aluminum levels, the accuracy 
obtained for 6 percent aluminum should be comparable to the first two cases. 

To illustrate some of the two-phase phenomena that occur, Figure 11 
shows a nozzle flow field for a different nozzle, the Space Shuttle SRB, which is 
basically a low’-altitude nozzle fired with a 16 percent aluminum solid propellant. 
At the exit of the nozzle the flow expands into a plume. A 10 pm particle, the 
innermost particle, does not turn with the flow; however, a smaller particle 
(3 pm) does turn some in the plume. The smallest particle shown is a 1.5 pm 
diameter particle. It initially shows an expansion and does tend to follow the 
plume boundary; however, the plume loses density and the ability to turn the 
particle back. At some point downstream, the small particle actually flows 
through the gaseous plume boundary. Thus, although calculations may show 
the plume boundary will not impinge upon a vent or orifice of the ship, particles 
passing through the plume boundary may actually be impinging or be captured by 
these ship orifices. 

3. Plume Afterburning Losses. The plume afterburning begins at the 
nozzle exit plane and is a function of the velocity or Mach number of the ambifpt 
air (Fig. 12). The air is entrained and finite chemical reactions occur; at some 
point the mixing layer reaches the nozzle axis, and then fully developed flow 
occurs downstream of that point. The species, the temperature, and the pressure 
vary as a function of the amount of mixing. One of the primary uses of this 
analysis is to determine how much air is entrained into the flow and to determine 
the dilution factor of the species that occur at the nozzle exit. This, in turn, 
allows one to determine the kinds of changes in species that occur due to the 
chemical reactions with the entrained air. 


C. Results for the HARPOON Booster Analysis 

The HARPOON booster nozzle was analyzed utilizing the two-dimensional, 
two-phase nozzle analysis program, TDDD [15] . The program was originally 
written by TRW Systems personnel and has, over the years, been modified by 
Science Applications, Inc. (SAI) personnel for NASA. The results of the nozzle 
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Figure 11. Nozzle-exhaust plume flow fields for a gas-particle 

system for a P P - 40S. 
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Figure 12. Plume afterburning schematic. 



analysis, in terms of static and total pressure at the nozzle exit, were used as 
Input parameters to the equilibrium program. The equilibrium program deter- 
mined the species present and their concentration (Fig, 12). A range of con- 
centration for eac ti species Is shown because the flow in the nozzle is neither 
uniform nor insentropic because of the two-phase flow. Every point in the flow 
field has a stagnation pressure which is dependent on the flow history of the 
particles in the neighborhood of the flow-field point. Due to the nonuniform 
flow in the nozzle, the static pressure varies radially across the exit plane. 

These two effects cause a species concentration variation across the exit plane. 
Table If gires the species concentration variation across the nozzle exit plane 
for the species of interest. 

Flow from the nozzle mixes and reacts with the ambient air at a finite 
rata. The LAPP ( 16} computer program calculates the amount of air mixed and 
the species concentrations radially from the nozzle axis to the freestream as a 
function of distance from the nozzle exit plane. Table 13 lists the finite-rate 
chemical reactions which are considered by the program for the HARPOON 
booster exhaust. Table 14 lists the results of the afterburning and mixing 
analysis. The species concentrations at the centerline of the nozzle and at the 
edge of the freestream are given for an axial station 120 ft downstream of the 
nozzle exit. This axial station was chosen because the temperature and velocity 
at the two racial stations have similar ..dues. The maximum concentration 
error due to thermodynamic data error is estimated for each species. Figure 13 
shows the static temperature of die centerline and edge as a function of distance 
from the nozzle exit. The edge static temperature rises more rapidly and to a 
higher value than the centerline. This is to be expected since the edge mixes 
more rapidly and sooner with the freestream. Figure 14 shows the dilution ratio 
or the relative amount of air that is entrained by the exhaust as a function of 
distance from the nozzle exit. Figure 15 shows the plume radial expansion 
versus distance. General conclusions that can be drawn are that because of the 
high nozzle exit pressure, 34.73 psia, the plume mixes rapidly with the free- 
stream even at a low vehicle velocity, afterburns near the nozzle exit, and 
expands and entrains approximately 500 times as much air as the exaaust 
effluent before the plume dissipates into the ambient atmosphere. 

Figure 16 is a photograph of an actual HARPOON missile flight. The 
photograph has been labeled to indicate some of the important points in the plume. 
For this particular nozzle, afterburning tai.t s p.ace very close to the nozzle exit. 
Fully developed flow (i.e. , the mixing region between the exhaust and the ambient 
air extends to the centerline) occurs approximately 20 ft aft of the nozzle exit. 

Air entrainment is essentially complete and the plume properties are uniform 
approximately 120 ft aft of the nozzle exit, A point of note is that every motor 
is different. The analysis performed for the HARPOON is not applicable to 
other motors; only the analysis technique is transfer ruble. 



TABLE 12. SUMMARY OF HARPOON BOOSTER EXHAUST 
SPECIES CONCENTRATIONS (ppm) 


© 

© 

© 

to 





CM 

§ 

CM 

CM 





• 

• 

oo 





© 

© 


© 





1 

| 

© 


- 


I 

r-1 l' - 


o 

rr 

■ 


\ - 

o ^ 

O , 

o 

% 

1 

o 

eg 

tH 

} 

O 

JC 

© 

o 


r-t 

1 


CM 

o 

o © 

o 

o 


LO 

1 



o u 


^ ^ 


© 







O 

o 

T— 1 


X 


CM 




© 

1 


"T 


T-H 

© 

t- 

CM 

CO 


CM 

1.0 

CM 


M 

1 

CM 

lO 



— 


» 

co 

CM 

1 

00 


l 

1 

t.O 

© 

1 7 
© 

1 

l0 

© 

| 

o 

CM 


CM 


o 

’TO 

© 



• 

• 

♦ 

♦ 

9 

• 

© 

# 

© 

CO 

© 


r-H 

o 

O 

© 

r-t 

o 


9 u u 


CO 

u 

r— 1 

§ 

S 

CM 

o 

r— < 


O 

u 

o 

CM 

r— 1 

8 

Cs 

o 

< 


< 

u 

o 

u 

U 

o 









TABLE 12. (Concluded) 
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TABLE 13. AFTERBURNING ANALYSIS (REACTIONS BEING CONSIDERED) 
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M represents any third body in the reaction. 





TABLE 14. HARPOON BC l TER EXHAUST EFFLUENT 
CONCEN ..ATIONS 


Species Concentrations (ppm) 

X = 120 ft downstream of nozzle exit 

P = 1800 psia 
c 

0. 1 sec after first ir ssile motion 

V_ = 412 K T . = 408 K 

£ edge 

V = 170 ft/sec V . =134 ft/sec 

edge 


Species 

Concentration 

Maximum Error 
( Percent) 

Centerline 

Edge 

HCl 

2932 

2824 

1 

h 2 o 

8290 

7985 

1 

CO 

33 

32 

1 

COz 

3519 

3390 

1 

Cl 

193 

186 

5 

ai 2 o 3 (s) 

3. 72 x 10 e fig/m 3 

3. 61 x 10 6 Mg/m 3 

1 

°2 

205 030 

205 210 a 

1 

N 2 

779 480 

779 870 a 

1 

Cl 2 

0.011 

0.011 

3 

NO 

0.04 

0.04 

1 

so 3 

1.8 

1.8 

1 

Cr 2 O s (S) 

1. 3 x 10 4 jug/m 3 

1.3 x 10 4 Mg/m 3 

- 

CuCl(L) 

8.4 x 10 4 Mg/m 3 

8.4 x 10 4 Mg/m 3 

1.8 


a. High concentration is from entrained air. 
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Figure 13. HARPOON booster afterburning analysis 
static temperature versus distance. 



Figure 14. HARPOON booster afterburning analysis 
dilution ratio versus distance. 
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V. CONCLUSIONS 


An analytical assessment of the exhaust effluents from the HARPOON 
booster was performed to idem the species present and determine their con- 
centrations in the exhaust plume. An error analysis was conducted to determine 
the potential uncertainty in the assessment. 

There were 95 candidate species (Table 3) considered in this investiga- 
tion. From the two-dimensional, two-phase thermochemical nozzle analysis 
it was determined that the thermochemical environment would support the 
formation of only 2S species (Table 12) in significant concentrations (above 
10" 6 ppm) at the exit plane of the nozzle. Using the finite-rate thermochemical 
analysis to account for the flow from the nozzle as it mixes and reacts with the 
ambient air, it was determined that the thermochemical environment reduced 
the numbers of significant species to 13 (Table 14) at 120 ft downstream of the 
nozzle. 




An error analysis of these results shows a maximum uncertainty in the 
concentrations of less than 1 percent in all but three cases. The uncertainty 

s 5, 3, and 1. 8 percent, respectively. Of these 


are so far below toxic 


for Cl, Clg, and CuCl^ was 

three species, the concentrations of CF and CuCl , v 

2 (e) 

limits that they can be neglected. The maximum uncertainty in Cl is ±10 ppm; 
however, in terms of other potentially toxic species present in these exhaust 
effluents, this small uncertainty can be neglected. 


This analysis clearly shows that the only species of potential concern 
are HC1 and A1 2 0 3 until the exhaust effluents become diffused into the atmosphere. 
Since monitoring A1 2 0 3 would be complex and could not be done in real time, it 
is recommended that only the HC1 be monitored. A recently developed analytical 
instrument, the Geomet, uses a chemiluminescent technique for the detection of 
HC1 or total chloride ion and is used by the Environmental Protection Agency [ 17) . 
This instrument 'uld be ideal for HCl detection. Thus, there is a minimum of 
analysis problems for the only gaseous species of concern. 


The study has accomplished its two main purposes. First, it has pre- 
dicted the exhaust species that are present and their concentration. Second, it 
has identified the species to be monitored. It has also found that there is no 
secondary species that can be used for monitoring purposes nor is there a need 
for one. The primary effluent of interest is easily, accurately, and inexpensively 
determined. 
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The analysis technique is genei’al and can be used for any solid rocket 
motor. Each motor has a different chamber pressure history, nozzle configura- 
tion, and propellant formulation. The results presented in this report are valid 
only for the HARPOON booster motor; they cannot be used to draw conclusions 
for any other motor. If tne exhaust effluents are of concern to personnel for 
assessment, each current and future motor should have similar analyses 
performed on it. For future motors, it would be most cost effective if the 
analyses were performed at an earlier stage of development and were made a 
permanent part of the standard development plan for a motor. 
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